clear all
d=readtable('graphdata.csv');
figure(1)
pl = plot(d.lbquant,d.sij','-k','linewidth',1.5);
hold on
p2 = plot(d.lbquant,d.lb','--^','color','red','markersize',4);
hold on
p3 = plot(d.lbquant,d.ub','-*','color','blue','markersize',4);
hold on
legend([pl p2 p3],'Estimated Choice Probability','Lower Bound','Upper Bound')
hold off
xlabel('Quantiles of the Lower Bound')
ylabel('Choice Probability')
ax = gca;
outerpos = ax.OuterPosition;
ti = ax.TightInset; 
left = outerpos(1) + ti(1);
bottom = outerpos(2) + ti(2);
ax_width = outerpos(3) - ti(1) - ti(3);
ax_height = outerpos(4) - ti(2) - ti(4);
ax.Position = [left bottom ax_width ax_height];
fig = gcf;
fig.PaperPositionMode = 'auto'
fig_pos = fig.PaperPosition;
fig.PaperSize = [fig_pos(3) fig_pos(4)];
axis([min(d.lbquant)-0.25   max(d.lbquant)+0.25   0  0.5])